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Abstract 

Complete wetting of geometrically structured substrates by one-component fluids with long- 
ranged interactions is studied. We consider periodic arrays of rectangular or parabolic grooves 
and lattices of cylindrical or parabolic pits. We show that the midpoint interfacial heights within 
grooves and pits are related in the same way as for complete wedge and cone filling. For sufficiently 
deep cavities with vertical walls and small undersaturation, an effective planar scaling regime 
emerges. The scaling exponent is —1/3 in all cases studied, and only the amplitudes depend on 
the geometrical features. We find quantitative agreement with recent experimental data for such 
systems. 

PACS numbers: 68.08.Bc,68.08.-p,05.70.Np 
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FIG. 1: Top view of the substrate geometries, (a) Quadratic lattice of identical pits (cylinders 
or paraboloids). (b) Periodic array of grooves (rectangular or parabolic). All cavities have finite 
depths D. 

The growing interest in device miniaturization has led to the emergence of various exper- 
imental techniques of tailoring the geometrical and chemical topography of solid surfaces at 
mesoscopic scales . Such nano-patterning of the surfaces may result in drastic changes 
of their wetting characteristics, which is important for technologies such as micro-fluidics 

n . q 

P| or the fabrication of super-hydrophobic or super- hydrophilic surfaces [2]. Experimental 
studies of complete wetting on sculptured surfaces jj, Q| demonstrate the strong influence 
of nanocavities on the adsorption behavior relative to that of flat substrates. Theoretical 
studies of adsorption in infinitely deep generalized wedges Q| predict geometry dependent 
wetting exponents. 

Although it is known that a non-planar topography of a substrate modifies its wetting by 
a fluid, recent studies have revealed surprising hidden symmetries, or so-called covariances, 
which relate various local adsorption properties for different substrate geometries. These 
covariances imply that different confining potentials can lead to identical local interfacial 
properties once external fields are suitably rescaled. For instance, for both long- and short- 
ranged forces, complete wetting at the apex of a wedge can be mapped onto critical wetting 
of a planar substrate with the apex angle playing the role of the contact angle Q , which maps 
one-to-one to temperature. In two-dimensional systems, critical wedge filling can be related 
to the strong-fluctuation regime of critical planar wetting |a|. Recently, a new example 
of geometrical covariance relating wedge and cone complete filling has been reported [^] 
showing that the equilibrium midpoint interfacial heights in a cone and a wedge obey 
the relation lj?\A[i, a) = /^(A/i/2, a) with a being the substrate tilt angle and A/z > 
the chemical potential deviation from liquid-vapor coexistence. This relation is valid for the 
leading behaviors of I in the limit A/i — > + . 
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Based on an effective interface Hamiltonian in the following we demonstrate that com- 
plete wetting of substrates patterned by periodic arrays of grooves or quadratic lattices of 
pits (see Fig. P), both of depths D, exhibits a geometrical covariance similar to the one de- 
scribed in Ref . |9| . We consider rectangular or parabolic grooves and cylindrical or parabolic 
pits, taking into account the long range of the intermolecular potentials. We observe four 
different scaling regimes: filling, post-filling, effective planar, and planar, with the neigh- 
boring regimes being separated by, as we call, A/i^ > A/i^ > A/x T . The aforementioned 
covariance relates the behavior of the midpoint wetting film thicknesses for all geometries 
and holds within an undersaturation range A//£ < A/i < A/Af(i?) (for A/j® see below, 
the superscripts p, g refer to pits and grooves, respectively) which we call the post-filling 
scaling regime. In the case of cylindrical pits or rectangular grooves, for A/i \ ApFjf^R) 
and for sufficiently large D/R, the analogue of capillary condensation occurs such that the 
pits or grooves are rapidly, but continuously, filled by the liquid. However, in the case of the 
parabolically shaped pits or grooves, A/j^f marks the crossover from the continuous power- 
law filling regime at A/i > A/x^f, to the post-filling scaling regime. We find numerically for 
R/a > 50 (see Fig. HJ, where a is a molecular length scale, Ap^f t (R) ~ R^ 1 ^ 5 with a small 
positive effective exponent 5, and Afjj U = 2Ap,^ u . If we denote the equilibrium interface 
height at the position of the symmetry axes of the cylindrical or parabolic pits as l p °\ and 
in the middle of the rectangular or parabolic grooves as l g , we obtain in the post-filling 
regime 



The scaling functions A P9 (x) do not depend on D and P, i.e., in this regime the midpoint 
interfacial heights increase upon decreasing undersaturation in the same way for an isolated 
cavity as for arrays of them; Ef is a molecular energy scale. 

For A/i < A/i^. the cavities are completely filled by the liquid and the equilibrium interface 
heights l Pt g(x,y) = l p>g become de facto independent of the lateral coordinates x = (x,y). 
In the case of the rectangular grooves or cylindrical pits, A/i^ marks the crossover from the 
post-filling scaling regime to the effective planar scaling regime within which the wetting 
behavior of geometrically patterned substrates can be mapped onto that of layered flat 
solids. The upper layer of those solids has a thickness D, and its composition is related to 



lj®(A»,R,P,D) = RA p>g ((An/s f )(R/a) 1+5 ) 




(1) 
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FIG. 2: Interfacial heights lp >g (measured from the top substrate layer z = 0) at the middle of the 
cylindrical and parabolic pits, l p ^ (dashed lines), and of the rectangular and parabolic grooves, l g ^ 
(solid lines), as a function of Afj, with C p = 1 and C g = 2. Rescaling the variables according to 
Eq.Q leads to data collapse. The inset shows the same data unsealed. The depth of the cavities 
is D = 500; pp(g) denotes parabolic pits (grooves). All length are measured in units of a. 

the geometrical parameters R and P. This results in the following scaling relations for the 
effective planar scaling regime: 

l p , g {An, R, P, D) = ($ Pi9 )^(A/i), (2) 

where l n is the thickness of the adsorbed liquid film on the corresponding non-patterned 
planar substrate. $ p = 1 — tt(R/P) 2 and $ 9 = 1 — 2R/P are the areal fractions of solid 
in the top layer, z — 0, of substrates with cylindrical pit and rectangular groove patterns, 
respectively. This result is reminiscent of the Cassie equation describing the apparent 
contact angle on chemically structured substrates. For D — > the width of the effective 
planar regime, i.e., the range of applicability of Eq. |2] vanishes. 

Finally, at A/i = Afi^ ~ with A{i n < the systems cross over to the planar 
scaling regime, in which the geometrical patterns are irrelevant. This crossover occurs only 
for long-ranged dispersion forces; we consider only them throughout. (For short-ranged 
interactions, instead, the growth of the wetting film would remain determined by the areal 
fraction of solid at z = for all film thicknesses.) In the case of the parabolic pits and 
grooves, we do not observe the effective planar scaling regime. There is rather an extended 
crossover region from the post-filling scaling regime to the planar one. 
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We base our calculations on the effective interface Hamiltonian description |h| 

H[l] = J J d 2 x (a lg y/l + (V/) 2 + ApApl + W(x, 0) , (3) 

where Z(x) denotes the local interfacial height which is measured from the plane z = where 
the substrate ends, oi g is the surface tension of the free liquid- vapor interface, Ap = p\—p g is 
the difference in number densities of the coexisting bulk phases, and W(x, I) is the effective 
interface potential. The functional in Eq. (J3J) can be derived systematically from microscopic 
density functional theory, which allows one to determine the explicit functional form of 
W(x, I) for a given substrate shape and its dependence on the jluid-fluid and substrate- 
fluid interaction potentials. We approximate the attractive parts of the pair potentials by 
0/,s( r ) = — ' £ } '(o" 2 ! + r 2 )~ 3 . The amplitude is chosen such that the integrated strength 
of <j)f equals that of the attractive contribution of the Lennard- Jones potential obtained by 
a strict application of the Weeks- Chandler- Andersen procedure. This leads to the effective 
interface potential 

W n (x,Z) = Ax/ n (x,Z), 
J n (x,Z) = [°°dz [ rfV(a 2 + |r-rf)" 3 , (!) 



where A = — 12 g^ 2 a 6 Ap(piEf — p s s s ) > is an effective Hamaker constant with p s as the 
number density of the substrate, and Q denotes the domain occupied by substrate particles. 
The effective interface potential of the planar substrate is W^{1 a) ~ An/(12l 2 ) so that 
h(Ap -»■ 0) = (A 7 r/(12ApA/i)) 1 / 3 . 

In the following we minimize the functional in Eq. Q numerically, which yields the 



equilibrium interface height I 



x) within mean-field theory which is valid in d = 3 for com- 



plete wetting ll| and filling ^| for the dispersion forces considered here. For all substrate 
geometries the midpoint heights l^ g (Ap) exhibit four different regimes. The first regime 
corresponds to the filling of the cavities. For the case of rectangular grooves and cylindrical 
pits, and for D/R large enough, a quasi abrupt, but still continuous filling of the cavities 
takes place at Ap = ApP^, which is shown in the inset of Fig. |21 A similar behavior has 
been reported earlier for an isolated rectangular groove |12]. We find that the locus Apfj^ 
of the filling transformation scales with the lateral size of the cavities as i? -1-5 , with an 
effective exponent 5 ~ 0.035 for both the rectangular grooves and the cylindrical pits. For 
the case of parabolic pits and grooves, complete filling of the structures is described by an 




FIG. 3: Wetting film thickness l g above grooves for small A/i. The symbols represent the midpoint 
interfacial height l g °\ The solid lines are obtained by minimizing Eq.© assuming /(x) = l g . On 
the present scales they are indistinguishable from those for film thicknesses /£(A/i) on layered and 
flat ersatz substrates, with a top layer of height D and an effective Hamaker constant A e = A& g , 
where A is the Hamaker constant of the solid without grooves. In the inset the vertical axis is 
rescaled according to Eq. @ leading to data collapse within an intermediate regime, the width of 
which increases as A/i£ — const x D~ 3 for D — > oo. The groove depth is D = 500. 

effective power law, /p°](A/i, R, D) ~ Afi~^ R ' D ^ valid for A/i > A/i?^. We find the values of 
the effective exponent 7 ranging from ca. 3.1, for R = 245cr, to ca. 2.0, for R = 50a, at a 
cavity depth D = 500cr. Moreover, the complete filling of the parabolic cavities obeys the 
covariance relation /p ^(A/i, R, D) = Z^(A/i/2, R, D). 

In the second, post-filling regime, i.e., for A/i G (A/i£, A/i^), the midpoint height for 
all patterns shows an almost linear dependence on A/i on normal scales. The morphologies 
of the liquid films still reflect the geometrical patterns, i.e., there are considerable lateral 
variations of the interfacial heights. We find that the slopes of the lp) g curves scale as R a with 
a « 2. Thus, combining this fact with the scaling of the filling chemical potential A/i^, we 
propose for the functions lp} g in the post-filling regime the scaling forms given by Eq. (Jljl. 
The scaling functions A Pt9 and the corresponding data collapse upon suitably rescaling the 
chemical potential are shown in Fig. |21 The scaling functions A p g in the post-filling regime 
do not depend on the cavity depth D and the pattern periodicity P. In the case of a single 
cavity, i.e., in the limit P — > 00, we obtain the same curves /p°g(A/i) as those presented in 
Fig. 121 

Below a certain value A/4 the interface height becomes de facto flat, denoted as l Pi9 
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for the respective patterns. In this case the functional in Eq. (jSJ) reduces to a function of 
minimum of which determines the thicknesses of the wetting films in this regime. 
Figure El shows the functions l g (Afi) for several values of R and P as well as the planar 
film thickness l n (Afi). The rectangular groove (and cylindrical pit, not shown) geometries 
lead to the same power law behavior as for a flat surface, Z Pi5 ~ A/i~ 1//3 , but with different 
amplitudes reflecting different effective Hamaker constants, which depend on the geometry 
of the patterns. We call this behavior the effective planar scaling regime. l g (Afi) for the 
parabolic grooves does not reveal the effective planar scaling regime, but gradually crosses 
over to the planar one, In order to find the geometrical dependence of the amplitudes of 
the scaling laws given above, we mimic sculptured substrates by layered and flat ersatz solids, 
with the effective interface potential W(l > a) « ^(A e /l 2 + (A-A e )/ (l+D) 2 ) Q. The first 
term is the effective interface potential of a flat semi-infinite solid with Hamaker constant 
A e and the second term is the correction due to the actual bottom part of the substrate, 
z G (—00, —D), with Hamaker constant A. For o <C I < D, to a good approximation one 
may ignore the bottom part of the substrate (i.e., the second term in W(l ^> a)) so that 
the amplitude is determined by the Hamaker constant A e . Consider a lateral unit cell ujq of 
this sculptured part of the substrate (z 6 [— D,0]), with u s as its domain occupied by the 
solid. By requiring W UWs (x, I) = {x, I), we obtain for the effective Hamaker constant 
A e = AI UoJs (-x,l)/I UoJQ (-x,l), where is the effective interface potential generated by 

the laterally homogeneous top layer of the layered and flat ersatz substrate; the symbol U 
denotes the union of domains. For I ^> cr the integrals I UuJi {i = s, 0) can be approximated 
as: 

f°° , f d 3 r> ( a 2 + llx-x'lPX 

In Eq. (JHJ) the leading term as function of I renders A e w AS s /Sq, where Si is the surface 
area of the domain uji D {(x, z = 0)}. Thus, we obtain A e g m A§ g for the rectangular grooves 
and A e p ~ A& p for cylindrical pits (see Eq. (J2J)). The thicknesses l^(Afi) of the wetting 
films on such layered ersatz substrates, which effectively correspond to arrays of grooves, are 
almost indistinguishable from the corresponding ones for l g . Therefore, we conclude that, 
for sufficiently thick wetting films, l g obeys the scaling relation given in Eq. (J2J). Indeed, 
after rescaling by the geometry-dependent factors (<3> s )~ 1//3 , the curves for l g collapse and, 
within the numerical precision, coincide with the curve for l n . This is shown in the inset 
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FIG. 4: Liquid adsorption at parabolic pits. Inset: T t (AT). Red line - theory; symbols and black 
lines - experimental data from Ref. J] (where T p is denoted as T c ). AT is proportional to Afi (see 
main text). 

of Fig. El In the case of a lattice of cylindrical pits we observe the same behavior of the 
corresponding l p curves (not shown). As l g reaches the value ~ D, i.e., at A/z = Afi n ~ D~ 3 , 
a crossover to the planar scaling regime takes place. In the planar scaling regime Afi < Afi n 
the geometrical patterns are irrelevant. As D decreases, Afi n eventually merges with A/i^., 
eliminating the effective planar scaling regime. 

Figure E] compares our results with the corresponding experimental data of Ref. 
in which the adsorption of methylcyclohexane (MCH) on a silicon substrate sculptured 
by a hexagonal lattice of parabolic pits has been studied with the following geometrical 
parameters: depth D m 200A; radius of the pits at the opening R « 123A; lattice constant 
P « 394A. As in Ref. U, we display the volume r = T p + T t of adsorbed liquid divided 
by the projected area and multiplied by the electron density of the bulk fluid as the sum of 
the amount T p adsorbed in the pit, and the amount T t adsorbed above the pit opening. For 
MCH we adopt the Lennard- Jones parameters a = 5. 511 A and Sf/ks = 446-fT 14|. Then the 
relation between the undersaturation Afi and the reservoir-substrate temperature difference 
(used in the experiment as a means to tune the deviation from liquid-vapor coexistence) is 
AT w 32.37^-fT. Fitting the film thickness on the planar substrate to the corresponding 
experimental curve fixes the value of the Hamaker constant. For the liquid- vapor surface 
tension of MCH we use a lg = 22.72^ for T = 30°C H- Adjusting the geometry of the 
patterns to the experimental ones we calculate T p and T t as functions of AT. We emphasize 
that at this stage there are no free parameters left in the model. The resulting curves 
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are shown in Fig. HJ For both quantities we obtain good quantitative agreement with the 
experimental results. In the filling regime (5K < AT < 8K) the calculated T P (AT) exhibits 
a power-law behavior T p oc AT _/3p with an effective exponent f3 p ~ 2.4. The apparent 
disagreement with the value (3 P ~ 0.76 used in Ref. jj] can be explained by the mislocation 
of the filling regime on the experimental T p (AT) curve due to the large error bars for the data 
points at large AT. The weak crossover at AT pa 8K corresponds to the disappearance of 
the interfacial meniscus, so that for larger values of AT the interface follows the shape of the 
substrate. We have found that for larger cavities this crossover becomes more pronounced 
and for sufficiently large D and R we obtain the exponent j3 p = 3.4. 

In summary, as function of undersaturation we have studied complete wetting of four 
classes of substrates structured by one- and two-dimensional periodic patterns of fixed depth 
and we have identified four scaling regimes. The filling and the post-filling evolutions of the 
interfacial profiles do not depend on the periodicity of the patterns, but are determined by 
a single isolated cavity. For sufficiently deep structures there exists a range of undersat- 
urations, in which the midpoint interfacial heights, obtained for different patterns, can be 
expressed in terms of a single scaling function. For small undersaturations, the single-cavity 
behavior crosses over to one dominated by the presence of many of them, for which the 
interfacial thickness increases as on a planar substrate, but characterized by an effective 
geometry-dependent Hamaker constant. Ultimately, for very small undersaturations the 
wetting film thickness becomes independent of the geometrical substrate structures. 
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